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摘 要 : 基于 Owens Valley Radio Observatory (OVRO) 40 米 望远镜 观测 数据 ， 我 们 收集 
T HE AA CGRaBS J0835+6835 在 15GHz 射 电波 段 的 数据 ( 约 12 年 ) 。 利 用 Lomb-Scargle 
periodogram (LSP)7 k 427042) RZA 4& (Weight Wavelet Z-transform, WWZ)Zr ik, 我 们 检测 
到 了 CGRaBS J0835+6835 存在 一 个 约 560 天 左右 的 准 周 期 震荡 (Quasi-periodic 
oscillation,QPO)。 我 们 运用 超大 质量 双 黑 洞 系 统 中 的 螺旋 喷 流 模型 ， 估 算 了 它 的 主 黑洞 质量 
M 的 范围 约 为 3. 7x10'Mo -3. 5X10 Mo。 为 了 探究 此 源 是 否 存在 聚 束 效应 ， 我 们 使 用 双 指 数 函 数 
拟 合 6 个 爆发 过 程 ， 发 现 光 变 多 普 勒 因子 8v 的 平均 值 约 为 10. 76. 这 表明 此 源 的 射电 辐射 存在 
明显 的 聚 束 效应 。 
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0 引言 


HEARS (Blazar) 是 活动 星系 核 (Active Galactic Nuclei, AGNs) 中 比较 特殊 的 一 个 
类 ， 是 最 活跃 的 天 体 之 一 。 其 辐射 主要 来 自 于 喷 流 中 的 高 能 粒子 ， 由 于 相对 论 聚 束 效应 ， 
使 得 其 可 以 在 伽 玛 波段 观测 到 辐射 活动 。 耀 变 体 可 分 平 谱 射电 耀 变 体 (FSRQ) 和 蝎 虎 座 BL 
型 天 体 (BL Lac) 两 个 子 类 "。 它 在 观测 上 的 特征 包括 高 光度 、 高 偏振 、 快 速 光 变 以 及 非 热 
辐射 等 。 能 谱 分 布 (Spectra Energy Distribution, SED) 时 现 出 双 峰 结构 ， 低 能 峰 是 
由 于 同步 辐射 产生 ， 高 能 峰 是 由 于 电子 的 逆 康 普 顿 辐射 过 程 或 者 强 子 辐射 产生 “”“"。 估 算 多 
普 勒 因子 的 方法 有 很 多 种 , 其 中 通过 射电 光 变 估算 多 普 勒 因子 (8v) 的 方法 是 比较 简单 的 一 种 
因而 被 广泛 应 用 "”” 。 光 变 研 究 能 探索 耀 变 体内 部 辐射 过 程 和 物理 机 制 ， 因 此 研究 天 体 的 光 
变 , 及 其 可 能 的 周期 性 ， 对 了 解 活动 星系 核 具 有 重要 意义 。 

ZB (ECGRaBS J0835+6835 (BWE 0830+6845) 是 一 颗 红 移 z=1. 414 的 天 体 (a=08 35 
47. 6080866785, 6=+68 35 11. 493386746) ”"”。Romani 等 人 测 到 这 个 源 在 R 波 段 的 流量 关 
0. 00013mJy”。Snellen 等 人 在 研究 GHz 峰 频 的 样本 源 时 , 给 出 了 此 源 在 325MHz、609MHz 和 5GHz 
等 不 同 频 率 下 的 流量 密度 ”。 本 文 利用 欧文 斯 谷 射电 天 文 台 (Owens Valley Radio Observatory, 
OVRO) 的 长 时 间 监 测 数 据 (15 GHz) 研究 CGRaBS J0835+6835 的 射电 光 变 ， 检 测 是 否 存在 准 周 
期 振荡 。 我 们 还 使 用 双 指 数 函 数 拟 合 爆发 过 程 ， 检验 是 否 存 在 明显 的 聚 束 效应 。 目 前 对 兆 变 
体 光 变 曲线 准 周期 振荡 进行 分 析 的 方法 主要 有 两 种 : Lomb-Scargle periodogram(LSP)"" 和 和 
加 权 小 波 Z2 变 换 (Weighted Wavelet Z Transform，W2T)""。 这 两 种 方法 非常 适合 处 理 非 等 
间隔 数据 。 本 文中 我 们 将 用 这 两 种 方法 对 次 变 体 CGRaBS J0835+6835 的 15GHz 射 电波 段 的 光 变 
周期 性 进行 研究 ， 并 对 得 到 的 结果 的 可 靠 性 和 置信 度 进行 评 佑 。 


* 基金 项 目 ; 国家 自然 科学 基金 (11863007) ; 云南 省 高 校 高 能 天 体 物理 重点 实验 室 资助 . 
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男 ， 
男 ， 


1 光 变 曲线 


如 图 1 是 来 自 于 Owens Valley Radio Observatory (OVRO) 40 米 望远镜 观测 的 耀 变 体 CGRaBS 
J0835+6835 在 15 GHz 射电 波段 的 光 变 曲线 ， 去 除 少 量 的 坏 点 后 ， 采 用 460 个 数据 点 ， 约 化 儒 
略 日 MJD 范 围 约 为 54908~58871。 我 们 采用 光 变 指数 A2 ”来 判断 天 体 的 活跃 程度 。 


A= J Soin (1) 
其 中 fs 和 fn 分 别 表示 流量 的 最 大 值 和 最 小 值 ， 光 变 指 数 A 值 越 大 表示 天 体 越 活 跃 。 通 过 
计算 可 以 得 到 CGRaBS J0835+6835 天 体 在 射电 波段 的 光 变 指数 为 0.98， 这 表明 CGRaBS 
J0835+6835 是 个 非常 活跃 的 天 体 。 
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图 1 蓝 线 表示 炮 变 体 CGRaBS J083546835 在 15 GHz 射电 波段 的 光 变 曲线 ， 红 线 表示 正弦 函数 对 光 变 曲线 的 大 致 拟 合 。 
Fig. 1 The blue line represents the light curve of the blazar CGRaBS J0835+6835 in the 15 GHz radio band, and the 


red line represents the rough fit of the sine function to the light curve. 


2 光 变 曲线 的 周期 分 析 
2.1 LSP 方 法 

检测 隐藏 在 噪声 中 的 周期 性 信号 是 天 文 数据 时 间 序 列 分 析 中 的 一 个 重要 目标 。 
Lomb-Scargle periodogram (LSP) 方法 可 以 对 非 均 匀 采 样 的 时 间 序 列 周期 图 进行 相位 修正 
处 理 ， 能 够 在 一 定 范 围 内 对 非 均匀 采样 的 时 间 间 隔 引起 的 误 报 周期 进行 修正 "… 。 因 此 ， 
Lomb-Scargle periodogram 方 法 能 很 好 的 寻找 隐藏 在 噪声 中 的 准 周期 振荡 光 变 。LSPF 方 法 的 
原理 是 基于 傅立叶 变换 ， 使 用 正弦 函数 模型 ， 基 本 公式 如 下 ": 

> (e) — x)cos olt, — r) »3 CC) 一 x)sin olt, 一 r) 


Px(@) = -二 十 22 
No 


S cos olt, — r) > sin? olt, - r) (2) 
j=1 Jsl 


T 
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其 中 + 为 对 应 时 间 t 的 相位 修正 ， 计 算 公 式 为 : 


xx sin?2w tj 
tan(2w7) = (3) 


2 
x cos*2wt; 


运用 (2) 式 有 几 个 优点 是 通常 的 离散 傅立叶 变换 所 不 具备 的 。 首 先 ， 公 式 中 包含 时 间 
项 t 使 得 周期 计算 不 随时 间 原 点 的 变化 而 变化 。 其 次 ， 这 种 形式 使 得 周期 图 分 析 完 全 等 价 于 


正弦 曲 
先 对 周 
Python 


示 准 周 


99%、99. 7% 置 信 度 。 图 3 中 的 绿 线 有 两 个 明显 的 峰值 分 别 为 562 天 和 1020 天 ， 并 且 两 个 峰 的 置 


线 对 数据 的 最 小 二 乘 拟 合 。 为 了 验证 LSP 方 法 计算 准 周期 的 正确 性 ， 如 图 2 所 示 ， 我们 
期 图 进行 寡 率 拟 合 ， 得 到 蒙特 卡 罗 方 法 计算 置信 度 的 所 需 参数 寡 率 指数 卫 9， 然 后 用 
编写 程序 计算 了 CGRaBS J0835+6835 天 体 的 置信 度 结果 ， 如 图 3 所 示 。 图 3 中 绿色 线 表 
期 图 ， 其 峰值 就 是 准 周期 结果 ; 蓝 色 、 红 色 、 黑 色 线 分 别 代 表 蒙 特 卡 罗 模 拟 的 95%、 


信 度 都 超过 了 99. 7%， 说 明 两 个 峰值 结果 都 是 可 靠 的 ， 但 是 1020 天 约 为 562 天 的 两 倍 ， 所 以 我 
们 取 最 小 正 周 期 即 562 天 作为 CGRaBS J0835+6835 天 体 的 准 周 期 结果 。 
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图 2 CGRaBS J0835+6835 在 射电 波段 的 对 数 寡 率 谱 。 


Fig. 2 CGRaBS J0835-6835 object logpower spectrum in radio band. 
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图 3 CGRaBS J083546835 在 射电 波段 的 QP0 的 置信 度 的 蒙特 卡 罗 模 拟 分 析 结 果 。 
Fig.3 The Monte Carlo simulation analysis results of the confidence level of the QPO of 


CGRaBS J0835-6835 in radio band. 


2.2 加权 小 波 Z 变 换 方 法 

小 波 分 析 (Wavelet Aanlysis) 是 一 种 新 的 分 析 方 法 ， 他 是 继 Fourier 分 析 之 后 纯粹 数 
学 和 应 用 数学 完美 结合 的 又 一 光辉 典范 。 利 用 Morlet 小 波 作为 小 波 母 函数 进行 变换 ， 复 
Morlet NIE Jg: 


w(t)= e "(go =e?) (4) 
Foster 定义 了 加 权 小 波 变换 (WWT) 5 : 
Wa CIV, S 
2V, 
其 中 
(sy Po) [Dot] — Xev€) [Ee] 
Nett = a Vy 二 一 E V. _ a a 
Pen E—e | m€eqU Eme | Ee 


由 于 小 波形 状 的 变化 ， 低 频 部 分 的 有 效 数 据 Nor 会 比 高 频 部 分 的 有 效 数 据 多 ， 因 此 
WWT 的 值 会 向 高 频 部 分 偏 移 ， 导 致 结果 出 现 侦 差 。 于 是 Foster 根据 他 所 提出 的 Z 统计 量 
定义 了 加 权 小 波 乙 变换 为 ! 

= (Nese y. (6) 
2(v. -v.) 


利用 (6) 式 可 以 得 到 如 图 4 的 准 周期 图 ， 左 图 是 时 间 周 期 平面 中 WZ 功率 彩色 缩放 的 分 
布 ， 右 图 是 时 间 平 均 功 率 ( 黑 色 曲 线 ) 作为 周期 的 函数 图 ， 其 中 蓝 色 虚线 表示 置信 概率 在 
99. 7% 和 95% 的 位 置 。 根 据 WZ 图 的 极 大 值 就 可 以 确定 光 变 曲线 的 周期 。 图 4 中 显示 了 一 个 
明显 的 极 大 值 560 天 且 置 信 概 率 超 过 了 99. 7%， 也 就 是 说 560 天 的 准 周 期 结果 是 可 靠 的 。 
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图 4 WWZT 方 法 对 CGRaBS J083546835 在 射电 波段 QP0 的 分 析 结果 。 
Fig.4 The WWZT test result for QPO search CGRaBS J0835+6835 at radio band. 
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3 多 普 勒 因子 分 析 


ZÉ WAT (8) 可 以 通 


两 个 量 都 是 难以 直接 观测 的 量 ， 因此， 有 多 月 
因子 的 方法 中 ， 通 过 射电 光 变 估算 多 普 勒 因子 (6v) 的 方法 相对 比较 简单 ， 因 而 被 广泛 应 用 。 
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过 喷 流 中 物质 流速 度 (V) 以 及 视角 (0) 两 个 本 征 参 数 定义 。 然 而 这 


日 间接 方法 估算 多 普 勒 因子 。 在 各 种 估算 多 普 勒 


+ 
asl Flare-1 0.25 Flare-2 
0.20} 9:29 
2 A 一 
全 ”015 上 Bor FIC 
g I Ht g ! 
c [4 
v v 
o p D 2 | 
x 010r x 0.10f} ~ 
E E | — 5 A | 
[n oe Sta uc 2^ t 
aa E. Ax 
905r. 0.05 + | im 
0.00 - asol 
—0.05 1 1 1 1 1 1 1 1 1 1 1 1 1 1| 
54800 54900 55000 55100 55200 55300 55400 55500 55600 55800 56000 56200 56400 56600 
MD MD 
Flare-3 | Flare-4 
0.20 - | 0.20 上 
) E 
> > 【 
> L = 2 n 
$ 0.15 ^ Poir vk 
E] L 5 / N 
5 \ 
E + 3 | } e 
5 N. 5 / N 
- m 
gi 010 n S V oior y v 
Fh / 
% , 
Ne z 机 
L B 7 * 
0.05 ae ` aosi * a 
e "i nl 
- E ss 
1 1 1 1 L L 1 1 1 1 1 1 1 1 
56400 56600 56800 57000 57200 57000 57100 57200 57300 57400 57500 57600 57700 57800 
MJD MJD 
,~ 
0.175 上  Flare-5 è K, o18-  Flare-6 
y $ 
0.150 上 ‘a * 0.16 - 
% 
/ + 
> / H > i 
2 0.125} / " i 0.14} fas 
E * E dis 
2 i S S + à " 
© 0.100} ^ S 0.12} yg \ 
x t t x "d \ + 
= 2 7 \ n 
z 4 T o10b } 5 M 
0.075 F r à "n P EY n 
4 - x: 
^ B t * " 
Fs ‘ 008- | | .--7 
0.050 H gr c Se A eee 
g^ S d + 
Pd a 0.06 - 
0.025} ^ 
1 1 
57600 57700 57800 57900 58000 58100 58200 58300 58200 58400 58600 58800 59000 
MD 


MD 


图 5 双 指 数 函 数 对 CGRaBS J0835+6835 的 6 个 爆发 (flare) 过 程 的 


NA. MID 为 约 化 儒 略 日 ， 单 位 是 天 。 


Fig.5 The double exponential function fits the curves of the 6 flares of CGRaBS J0835+6835. MJD is Modified Julian 


Day, the unit is day. 
文献 [19-20] 中 详细 介绍 了 亮 温度 和 双 指 数 函 数 公式 及 估算 多 普 勒 因子 的 过 程 。 双 指数 函数 
公式 如 下 : 


for tto 
F(t) = F.+Fo(e & teta) (7) 


我 们 利用 亮 温 度 公式 估算 CGRaBS J0835 


+6835 天 体 的 射电 亮 温度 约 分 布 在 3.7X10 ~2.3X 


10" K， 超 过 首 康 普 顿 极限 温度 或 平衡 亮 温度 几 个 量 级 ， 表 明 此 源 的 射电 辐射 存在 明显 的 聚 


束 效应 。 双 指数 函数 拟 合 的 6 个 爆发 过 程 如 图 
步 估 算 对 应 的 光 变 多 普 勒 因子 。 每 个 爆发 过 程 的 拟 合 参数 如 表 1， 第 1 列 为 
范围 ， 第 2 列 为 拟 合 所 得 的 约 化 最 小 残 差 平方 和 pf， 


约 化 儒 略 日 
第 3 列 为 归 


拟 合流 量 ， 


进 一 


5 所 示 。 类 似 文献 [21] 选 取 平衡 亮 温 度 作为 内 


详细 的 定义 见 文献 [22]; 


x +r Ps \y 、 i-fiV? -二 ` M Y ES v 
一 化 后 的 卡 方 值 ， 其 定义 为 zz/dof， 妇 = Yu. (PSH), Jtr Pi 为 实测 流量 ，f 为 


o 为 残 差 绝 对 值 的 标准 偏差 ， 


o 


自由 度 dof=N—5; 第 4 列 为 拟 合 峰 值 对 应 的 约 化 


儒 略 日 及 误差 ; 第 5 列 为 基底 流量 及 误差 ; 第 6 列 为 爆发 的 幅度 及 误差 ; 第 7 列 和 第 8 列 分 


别 为 指数 上 升 时 标 和 下 降 时 标 及 误差 ; 第 9 列 为 每 个 爆发 过 程 所 对 应 的 多 普 勒 因子 。5v 的 值 


分 布 在 8. 87+0. 68-14. 4540. 23， 平 均值 约 为 10. 76。 


表 1 爆发 过 程 的 拟 合 结果 
Tablel The results of flare fitting 


B xeldof 
peak range to Fc Fo tr ta dv 
(2) (3) 

(1) (4) (5) (6) (7) (8) (9) 
54908.9-55366.2 1.52E-04 1.54 55062.97+26.79 0.06+0.04 0.26+0.03  140.51449.16 69.25+24.09 9.63+0.34 
55623.2-56414.8 3.98E-04 2.94 55978.50+9.63 0.07+0.04 0.33+0.11 152.17+48.32 98.21+56.41 8.87+0.68 
56487.6-57146.8 3.14E-04 1.63 56909.21+27.48 0.14+0.02 0.31+0.17 146.58+20.94 64.17444.78 10.4340.35 
57186.7-57671.5 1.43E-04 1.95 57336.36526.94 0.24+0.06 0.33+0.03 95.97+28.51 123.84+24.31 9.36+0.12 
57717.3-58174.0 4.23E-04 1.15 57978.93+11.71 0.11+0.01 0.3640.03 123.38+43.37 53.92+19.23 11.77+0.84 
58257.7-58724.5 1.08E-04 2.79 58589.75+18.08 0.07+0.04 ”0.11+0.01 33.51+38.41 128.52+47.14  14.4540.23 


4 讨论 与 结论 


通过 收集 欧文 


周期 约 为 562 天 和 1020 天 ， 并 
故 可 以 认为 CGRaBS J0835+6835 的 准 周 


斯 谷 射电 天 文 台 的 长 时 间 监 测 数 据 ,我们 利 ) 
法 对 CGRaBS J0835+6835 的 射电 波段 光 变 


JLSP 方 法 和 加 权 小 波 Z 变 换 方 
| 线 的 准 周期 振荡 进行 分 析 。LSP 方 法 分 析 得 到 的 准 


560K, 


其 置信 和 度 都 大 于 36 (99.7%) ，1020 天 大 约 是 562 天 的 两 倍 ， 
期 约 为 562 天 。 加 权 小 波 7Z 变 换 方法 得 到 的 准 周 期 约 为 


置信 和 度 也 超过 了 3o (99. TX) 。 两 种 检测 准 周期 的 方法 得 到 的 结果 都 大 致 约 为 560 


天 ， 检 验 的 结果 相互 印证 ， 所 以 我 们 认为 CGRaBS J0835+6835 天 体 的 准 周期 约 为 560 天 。 同 


时 ， 我 们 利 


j 双 指数 函数 拟 合 了 6 个 爆发 过 程 ， 


得 到 5v 的 值 分 布 8. 87+0. 68-14. 45+0. 23， 平 


均值 约 为 10.76。 结 果 表 明 CGRaBS J0835+6835 的 射电 辐射 存在 明显 的 聚 束 效应 ， 支 持 相 对 论 


性 喷 流 模型 


已 有 很 多 物理 


。 光 变 体 长 周期 光 变 的 物理 模型 现在 仍然 是 悬而未决 的 问题 ,为 了 解决 这 些 问题 ， 


模型 被 提出 ， 例 如 : 双 黑 洞 模型 下、 螺旋 喷 流 模型 中 和 薄 吸 积 盘 理论 等。 
根据 CGRaBS J0835+6835 天 体 产 生 的 560 天 准 周 期 , 其 物理 过 程 可 能 是 喷 流 的 


螺旋 运动 ”。 


Rieger 在 2004 年 给 出 了 观测 准 周期 三 实际 物理 驱动 周期 已 之 间 的 关系 公式 为 : 


fb iA Ia 
C hinaA IVA 


2 
Yo 


Lež (8) 
其 中 红 移 z=1. 414”，y, 是 体 洛 伦 兹 因子 约 等 于 7.5””。 基 于 (7) 式 可 以 得 到 实际 物理 驱 
动 周期 P35. 75yr。 另 外 ， 如 果 y 志 15 VEU ZU, P2143. 00 Yr。 如 果 主 黑洞 和 次 黑洞 
之 间 的 质量 比 1/3， 称 之 为 “major merger” (EFE) 。 如果 3<R<10', 称 之 为 “minor 
merger”( 次 并 合 )”“”。 无 论 质量 比如 何 ， 主 黑洞 的 质量 都 可 以 用 如 下 公式 估计 : 

Mx P3 R^10M,, (9) 
其 中 PLEA RAW, IF SMBBH 系统 的 “ 主 并 合 ”， 质 量 比 可 以 假设 为 让 3/2。 把 我 们 的 
参数 代入 到 (8) 式 中 ,得 到 CGRaBS J0835+6835 的 主 黑洞 质量 大 约 为 M-3. 7xl0Mo。 同 理 ， 
如 果 使 用 y=15 作为 参数 ， 次 并 合 的 SMBBH 系统 的 主 黑 洞 质量 为 1 =3. 5x10 Me， 故 我 们 取 
CGRaBS J0835+6835 的 主 黑洞 质量 V 的 范围 约 为 3. 7x10 Me~3. 5x10°Mo» 


P, = 
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Radio quasi-period oscillation analysis of blazar 


CGRaBS J0835+6835 


Gong Yunlu, Yi Tingfeng", Yang Xing, Chang Xin’, Li Yong 
(1.School of physics and electronic information, Yunnan Normal University, Kunming 650500, China, yitingfeng98 € 163.com; 2, 


Guangxi Key Laboratory for the Relativistic Astrophysics, Nanning 530004, China) 


Abstract: Based on the Owens Valley Radio Observatory (OVRO) 40m telescope observation 
database, we collected data from the blazar CGRaBS J083546835 in the 15 GHz radio band 
(approximately 12 years). Using the Lomb-Scargle periodogram (LSP) method and the Weight 
Wavelet Z-transform (WWZ) method, we detected a quasi-period oscillation of about 560 days in 
the CGRaBS J0835+6835 object. We used a spiral jet model in the supermassive dual black hole 
system to estimate the mass range of its primary black hole M range is 3.7x10°Mo to 3.5x10" Me. 
To investigate whether there is a beam effect in this source, we used a double exponential function 
to fit six burst processes and found that the average value of the variability Doppler factor dv is 
about 10.76. This result indicates there is a significant beam effect in the radio radiation of this 
source. 

Keywords: Blazar; CGRaBS J0835-6835; LSP Method; WWZ Method 


